Distribution of dwell times of a ribosome: 
effects of infidelity, kinetic proofreading and ribosome crowding 
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Ribosome is a molecular machine that polymerizes a protein where the sequence of the amino acid 
residues, the monomers of the protein, is dictated by the sequence of codons (triplets of nucleotides) 
on a messenger RNA (mRNA) that serves as the template. The ribosome is a molecular motor that 
utilizes the template mRNA strand also as the track. Thus, in each step the ribosome moves forward 
by one codon and, simultaneously, elongates the protein by one amino acid. We present a theoretical 
model that captures most of the main steps in the mechano-chemical cycle of a ribosome. The 
stochastic movement of the ribosome consists of an alternating sequence of pause and translocation; 
the sum of the durations of a pause and the following translocation is the time of dwell of the 
ribosome at the corresponding codon. We derive the analytical expression for the distribution of the 
dwell times of a ribosome in our model. Whereever experimental data are available, our theoretical 
predictions are consistent with those results. We suggest appropriate experiments to test the new 
predictions of our model, particularly, the effects of the quality control mechanism of the ribosome 
\^ and that of their crowding on the mRNA track. 
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PACS numbers: 87. 16. ad, 87.16.Nn, 87.10.Mn 



I. INTRODUCTION 

The primary structure of a protein consists of a sequence of amino acid residues linked together by peptide bonds. 
^3 ■ Therefore, a protein is also referred to as a polypeptide which is essentially a linear hetero-polymer, the amino acid 
residues being the corresponding monomers. The sequence of the amino acid residues in a polypeptide is dictated by 
that of the codons, each of which is a triplet of nucleotides, on a messenger RNA (mRNA) that serves as the template 
for the polypeptide. The actual polymerization of each protein from the corresponding mRNA template is carried 
i-C , out by a macromolecular machine called ribosome [H-Q and the process is referred to as translation (of genetic code). 
The polymerization of protein takes places in three stages that are identified as initiation, elongation (of the protein) 
\ and termination. 

The ribosome also utilizes the template mRNA as a track for its own movement; it steps forward by one codon at 
a time while, simultaneously, it elongates the growing polypeptide by one amino acid monomer. Therefore, ribosome 
is also regarded as a motor that, like other molecular motors, takes input (free-) energy from the hydrolysis of 
a nucleotide tri-phosphate to move along a filamenous track |4j. In fact, a ribosome hydrolyzes two molecules of 
Guanosine tri-phosphate (GTP) to move forward by one codon. 

Enormous progress has been made in the last decade in the fundamental understanding of the structure, energetics 
and kinetics of ribosomes [5l-fl4j . Recent single molecule studies of ribosomes fl5l - [2l| has thrown light on its operational 
mechanism. 

In single molecule experiments, it has been observed that a ribosome steps forward in a stochastic manner; its 
stepping is characterized by an alternating sequence of pause and translocation. The sum of the durations of a pause 
and the following translocation defines the time of a dwell of the ribosome at the corresponding codon. The time of 
dwell of a ribosome varies randomly from one codon to another. This randomness arises from two different sources: (i) 
intrinsic fluctuations associated with the Brownian forces as well as the low of concentrations of the molecular species 
involved in the chemical reactions, and (ii) extrinsic fluctuations caused by the inhomogeneities of the sequence of 
nucleotides on the template mRNA [22J. Because of the sequence inhomogeneity of the mRNA templates used by 
Wen et al. (23[, the dwell time distribution (DTD) measured in their single- molecule experiment reflects a combined 
effect of the intrinsic and extrinsic fluctuations on the dwell time. In contrast, in this paper we focus on situations 
(to be explained in detail in the following sections) where the randomness of the dwell times arises exclusively from 
intrinsic fluctuations. 

The probability density f(t) of the dwell times of a ribosome, measured in single-molecule experiments [2j|, has 
been compared with the corresponding data obatined from computer simulations 24]. It has been claimed that 
the data do not fit a single exponential thereby indicating the existence of more than one rate-limiting step in the 
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mechano-chemical cycle of each ribosome. In fact, the best fit to the simulation data was achieved with five different 
rate-determining steps (24|. 

A systematic analytical derivation of the DTD of the ribosomes was presented recently [25[ from a kinetic theory 
of translation (2|| that also involves essentially five steps in the mechano-chemical cycle of a ribosome. However, 
the model developed in ref.|26| ignores some of the key features of the mechano-chemical cycle of an individual 
ribosome during the elongation stage. For example, a ribosome deploys an elaborate proofreading mechanism to 
select the correct amino acid dictated by the template (and to reject the incorrect ones) to ensure high translational 
fidelity. Nevertheless, occasionally, translational errors take place. In this paper we extend Basu and Chowdhury's 
original model [26] by capturing the processes of proofreading and allowing for the possibility of imperfect fidelity of 
translation. Moreover, the identification of the mechano-chemical states as well as the nature of the transitions among 
these states is revised in the light of the empirical facts established in the last couple of years. Using this revised and 
extended kinetic model of translation [27| . we analytically calculate the probability density f(t) of the dwell times of 
ribosomes. 

The DTD derived in this paper is qualitatively similar to that observed by Wen et al. in their single molecule 
experimen ts [ 23l| . However, because of the sequence inhomogeneity of the template mRNA used in the experiment of 
Wen et al.[23|. their data cannot be compared quantitatively with our analytical expression for f(t). Therefore, we 
propose a concrete experimental set up required for a quantitative testing of our theoretical predictions. However, 
till such experiments are actually carried out, the results reported here will continue to provide insight into the 
mechanistic origin of the qualitative features of the DTD. Moreover, attempts are being made to extend our model 
to capture sequence inhomogeneities of mRNAs and to calculate the corresponding DTDs of ribosomes analytically. 

Interestingly, the inverse of the mean dwell time is the average velocity V of the ribosome motor. Since, in our 
model, the ribosome is not allowed to back track on the mRNA template, V is also the mean rate of elongation of 
the polypeptide. We show that V satisfies an equation that resembles the Michaelis-Menten equation for the average 
rate of simple enzymatic reactions [28j . 

It is well known that most often a large number of ribosomes simultaneously move on the same mRNA track each 
polymerizing one copy of the same protein. This phenomenon is usually referred to as ribosome traffic because of 
its superficial similarity with vehicular traffic on highways [HI, [29l - |4]| In this paper we also report the effect of steric 
interactions of the ribosomes in ribosome traffic on their DTD. 



II. THE MODEL 

We develope here a theoretical model for the polymerization of a protein by ribosomes using an artificially synthe- 
sized mRNA template that consists of a homogeneous sequence (i.e., all the codons of which are identical). In the 
surrounding medium, two species of amino-acids are available only one of which is the correct one according to the 
genetic code. The ribosome deploys a quality control system that rejects the amino acid monomer if it is incorrect. If 
this quality control system never fails, perfect fidelity of translation would result in a homopolymer whose constituent 
amino acid monomers are all identical. However, occasionally, wrong amino acid monomer escapes the quality control 
system. Such translational errors, whereby a wrong amino acid monomer is incorporated in the elongating protein, 
results in a hetero-polymer. We'll study the effects of the quality control system on the DTD of the ribosomes. 

A ribosome consists of two interconnected subunits which are designated as "large" and "small" (see fig ID . The 
small subunit binds with the mRNA track and decodes the genetic message of the codon whereas the polymerization of 
the protein takes place in the large subunit. The operations of the two subunits are coordinated by a class of adapter 
molecules, called tRNA (see figO}. One end of a tRNA helps in the decoding process by matching its anti-codon with 
the codon on the mRNA while its other end carries an amino acid; in this form the complex is called an amino-acyl 
tRNA (aa-tRNA). 

The three main stages in each mechano-chemical cycle of a ribosome, shown in fig|21 are as follows: (i) selection 
of the cognate (i.e., correct) aa-tRNA, (ii) formation of the peptide bond between the amino acid brought in by the 
selected aa-tRNA and the elongating protein, (iii) translocation of the ribosome by one codon on its track. However, 
some of these steps consist of important sub-steps. Moreover, the aa-tRNA selected (erroneously) by the ribosome 
may not be the cognate tRNA; this leads to a branching of pathways. Such branching, in turn, gives rise to the 
possibility of more than one cyclic pathway for a ribosome in a given cycle. 

Let us begin with the state labelled by "1" ; both the E and A sites are empty while the site P is occupied as shown 
in figEl An incoming aa-tRNA molecule, bound to an elongation factor EF-Tu, occupies the site A; the resulting state 
is labelled by "2". The transition 1—^2 takes place at the rate tu a . However, not all incoming aa-tRNA molecules 
are automatically selected by the ribosome. In order to ensure translational fidelity [42Tj45j the ribosome deploys a 
quality control mechanism to ensure that the aa-tRNA selected is, indeed, cognate (i.e., carries the correct amino acid 
as dictated by the mRNA template). A non-cognate tRNA is rejected on the basis of codon- ant ico don mismatch; the 
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FIG. 1: Cartoons representing the main molecules and machinery involved in the process of translation. 




FIG. 2: Pictoral depiction of the three main stages in the chemo-mechanical cycle of a single ribosome (see the text for details). 

corresponding transition 2 — > 1 takes place at the rate uj t \. However, if the aa-tRNA is not a non-cognate one, the 
GTP molecule bound to it is then hydrolyzed to GDP causing the irreversible transition 2 — > 3 at the rate cj/a- At 
this stage, kinetic proofreading can reject the tRNA if it is a near-cognate one; the transition 3 — > 1 captures this 
phenomenon and takes place at the rate ui r 2- In our model we do not explicitly treat those sub-steps in which EF-TU 
and the products of the hydrolysis of GTP leave the ribosome. We'll later utilize the fact that uj a is proportional to 
the concentration of the aa-tRNA molecules in the surrounding medium. 

If the selected aa-tRNA is not rejected, it does not necessarily imply that it is cognate because, occasionally, even 
non-cognate or near-cognate aa-tRNA escapes the ribosomes quality control mechanism. The amino acid supplied by 
the selected aa-tRNA is then linked to the growing polypeptide by a peptide bond; this peptidyl transferase activity 
of the ribosome elongates the polypeptide by one monomer. The correct amino acid is incorporated at a rate w p 
whereas a wrong amino acid is incorporated at a rate fl p thereby giving rise to the two different branches of the 
mechano-chemical cycle of the ribosome; the third branch 3 — ► 1, however, ends a "futile" cycle. The only difference 
between the states 4 and 4* is that the last amino acid in the polypeptide is correct in 4 but incorrect in 4* . While the 
polypeptide gets elongated by one amino acid, a fresh molecule of GTP enters bound with an elongation factor EF-G. 
The arrival of GTP-bound EF-G is not treated explicitly in our model. Alternatively, lo p (and f2 p ) is an effective rate 



FIG. 3: Pictoral depiction of the main steps in the chemo- mechanical cycle of a single ribosome (see the text for details). 



constant that accounts for both the polypeptide elongation and the arrival of the GTP-bound EF-G. 

Next, spontaneous Brownian (relative) rotation of the two subunits of the ribosome coincides with the back-and- 
forth transition between the so-called "classical" and "hybrid" configurations of the two tRNA molecules. [4614561 ] In 
the classical configuration, both ends of the two tRNA molecules correspond to the locations of P and A sites. In 
contrast, in the hybrid configuration the ends of tRNA molecules interacting with the large subunit are found at the 
locations of E and P sites, respectively. The forward transition 4 — s> 5 takes place at the rate w&/ along the correct 
branch (and at the rate tlbf along the wrong branch) whereas the reverse transition 5 — > 4 takes place at the rate 
LJbr along the correct branch (and ttbr along the wrong branch). The reversible transition 4^5, which is caused by 
spontaneous Brownian fluctuations, does not need any free energy input from GTP hydrolysis. 

Finally the hydrolysis of GTP drives the irreversible transition 5 — > 1 which involves the translocation of the 
ribosome on its track by one codon and, simultaneously, that of the tRNAs inside the ribosome by one binding site 
following which the deacetylated (i.e., denuded of amino acid) tRNA exits from the E site. The rate of this transition 
is LUh2 or depending on whether the system is following the correct branch or the wrong branch. In our model, 
the transition 5 — > 1 captures the hydrolysis of GTP, departure of the products of hydrolysis alongwith EF-G as well 
as the final exist of the deacetylated tRNA molecule from the E site. 

The rate constant ui a can be made identical for both species of amino acid monomers by maintaining their concen- 
trations in the medium appropriately. The rate constant uj r \ refers only to the wrong amino acids (i.e., non-cognate 
tRNA) because, we assume, cognate tRNA is not rejected at all. Since the transition 2 — > 3 accounts for the hydrolysis 
of a GTP molecule by the GTPase EF-Tu, the corresponding rate constant is u>hi, irrespective of the identity of the 
attached aa-tRNA. In kinetic proofreading, the rate of rejection of the non-cognate tRNA is much higher than that 
of a cognate tRNA. For the sake of simplicity, we assume that the cognate tRNA is not rejected at all. Therefore, 
the rate constant uj r2 referers to the rejection of only non-cognate tRNA. For the remaining steps of the mechano- 
chemical cycle, the rate constants are w p , wt,y, Ub r and o;^ 2 provided a cognate tRNA has been selected finally. On the 
other hand, the corresponding rate constants are Q p , Qbf, &br and flh2, respectively, if a non-cognate tRNA escapes 
rejection by the quality control system of the ribosome. Since the last step involves not only hydrolysis of GTP by 
the GTPase EF-G, but also translocation of the ribosome and the tRNA molecules, the rate constants ujya and VLki 
need not be equal, in general. 

Thus our model is an extension of the generic models for molecular motors based on stochastic chemical kinetic 
approach which was pioneered by Fisher and Kolomeisky [53, HH . Following Hill [5^| the mechano-chemical kinetics 
of a ribosome in our model can be regarded as a composite of the four cycles shown in FigUJ 

All the numerical data generated from the analytical expressions for the graphical plots have the following set 
of values of the rate constants (except in the figures where fit) has been plotted for several different values of the 
parameters ui p and Lo r2 ) : cu a — 25s" 1 , u rl = 10s _1 , u>hi = 25s -1 , O p = 40s _1 , u>bf = ^br = 25s" 1 , ftbf = ^hr = 10s _1 , 
ujh2 = 25s -1 , Qh2 = 10s -1 . The value of uj a = 25s _1 is identical to that used in our earlier papers (see [25|,[26[ and the 
references therein). For the purpose of plotting our results, the magnitude of the other parameters have been chosen 
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FIG. 4: The full mechano-chemical kinetic of a ribosome during the elongation state, as modelled in Fig(3] can be regarded 
as a composite of at least four different cycles shown in (a)-(d). Among these, both (a) and (b) are unproductive in the sense 
that these lead to neither elongation of the polypeptide nor forward stepping of the ribosome. However, (b) is "futile" because 
it dissipates the free energy input from the hydrolysis of a GTP molecule whereas no fuel is wasted in (a) . In contrast, both 
(c) and (d) lead to elongation of the polypeptide and forward stepping of the ribosome by hydrolysing two molecules of GTP; 
however, (c) incorporates the correct amino acid whereas (d) incorporated a wrong amino acid into the growing polypeptide. 

to be comparable to that of Lo a . The values of some of the parameters, and the ranges over which some of these have 
been varied, may be unrealistically high or low. But, this deliberate choice has been motivated by our intention to 
demonstrate graphically the interplay of various kinetic processes in translation. 



III. DWELL TIME DISTRIBUTION 



For the convenience of mathematical calculations, we assume that, after reaching the state 5 (or 5*) at location j 
the system makes a transition to a hypothetical state 1 at location j + 1 which, then, relaxes to the state 1, at the 
same location, at the rate 8. At the end of the calculation our model is recovered by setting 6 — > oo. 

Suppose Pfj,(j, t) denote the probability at time t that the ribosome is in the "chemical" state /j, and is decoding the 
j-th codon. Let us use the symbol P\{j + to denote the probability of finding the ribosome in the hypothetical 
state 1 at time t while decoding the j + 1-th codon. The time taken by the ribosome to reach the state 1 at j + 1, 
starting from the initial state 1 at j, defines its time of dwell at the j-th codon. Since, in this context, all the 
"chemical" states, except 1 refer to the J-th codon while 1 corresponds to the j + 1-th codon, from now onwards we 
drop the site index j (and j + I) to keep the notations simple. 

The master equations governing the time evolution of the probabilities P^(t) can be written as 



dPi{t) = 
dt 

dP 2 (t) 
dt 

dPsjt) _ 
dt 

dPi(t) 
dt ' 



-w a P x (t) + oj rl P 2 (t) + u r2 P 3 (t) (1) 

= UJaP.it) - [un + 0J hl )P 2 (t) (2) 

Lo hl P 2 (t) - (oj p + n p + uj r2 )P 3 (t) (3) 

= LJ p P 3 (t) - u bf Pi{t) + uj bT P 5 (t) (4) 
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FIG. 5: The probability density of the dwell times are plotted for a few different values of the parameter uj p . 

dP 5 {t) 



dt 

dPjjt) 
dt 

dpqt) 

dt 



u b fP 4 {t) - {w m + u br )P 5 (t) 

= n p p 3 (t) - n bf p*(t) + n br p;(t) 
n bf pz(t)-(n h2 + n br )p r ;(t) 



dPijt) 
dt 



tu h2 p 5 {t) + n h2 p;(t) 



Because of the normalization condition 



p^t) + put) + p B *(t) + Pt(t) = i 

not all the equations CQ)-© above are independent of each other. 

For the calculation of the dwell time, we impose the initial conditions 

Pi(0) = 1, and P 2 (0) = P 3 (0) = P 4 (0) = P 5 (0) = P 4 *(0) = P 5 *(0) = P(0) = 



(5) 
(6) 
(7) 

(8) 
(9) 



(10) 



Suppose, f(t) is the probability density of the dwell times. Then, the probability of adding one amino acid to the 
growing polypeptide in the time interval between t and t + At is f(t)At where 



fit) - = ^ 2 p 5 (i) + n h2 p;(t) 

The calculation of the DTD [6fjl - l65l | is essentially that of a distribution of first-passage times [if 



(11) 
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FIG. 6: The probability density of the dwell times are plotted for a few different values of the parameter u} r 2- 




FIG. 7: The probability density of the dwell times are plotted for a few different values of the coverage density p cov of the 
ribosomes. 
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The exact probability density of the dwell times is given by 
/(<) = 
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where cj! ,cj 2 and w 3 are solution of the cubic equation 

W 3 — UJ 2 (iO r l + Uh\ +td a + W r2 + UJ p + ftp) + Uj(uJhl^a + UJ r2 UJrl + UJ r 2^hl + W r2 CJ a 
+ UJpUJrl + UJpUJhl + UpLO a + flpU! r l + VtpbJhl + £^pW a ) — VlpLOhlUJa + LOpU!hlU! a = 0, 

w 4 and W5 are the solution of the quadratic equation 

- Uj(uj h2 + LUbr + U b f) + Ld h2 UJ b f = 

and f2 4 and fis are the solution of the quadratic equation 

n 2 - Q{n h2 + n br + n bf ) + n h2 n bf = o 



(12) 

(13) 
(14) 
(15) 



Some of the details of this derivation are given in the appendix. Note that the problem of determining the rates 
uji (i — 1,2,3,4,5) in terms of the rate constants of the model is similar to that of expressing the normal modes of 
vibration of a set of coupled harmonic oscillators. It is the "backward" reactions in the mechano-chemical cycle of 
our model which play the role of coupling of the harmonic oscillators. For example, if Q p — and simultaneously, 
u! r i = oj r2 — uj b r — 0, the expressions for (i = 1, 2, 3, 4, 5) reduce to the simple form w± = uj a , u 2 — 0JhU w 3 = w p 5 
w 4 = u)bf, and uj 5 = uj h2 . 

The distribution (fT2|) is plotted in Figl5]for a few different values of the parameter ui p . Note that 



4> = UJp/(uJ p + tt p ) 



(16) 



is a measure of the fidelity of translation. Therefore, increasing ui p , keeping fl p fixed, enhances translational fidelity. 
Moreover, a higher uj p also corresponds to faster peptidyl transferase reaction. Therefore, the trend of variation of 
the most probable dwell time with uj p is consistent with the intuitive expectation that the slower is the peptidyl 
transferase reaction, the longer is the dwell time. 

The kinetic parameter w r2 is a measure of the rate of rejection of the aa-tRNA by kinetic proofreading. Consequently, 
the effect of the variation of the uj r2 on the most probable dwell time is opposite to that of uj p (see FigEJ; the higher 
is the frequency of tRNA rejection by kinetic proofreading, the longer is the dwell time. 



IV. MEAN RATE OF POLYMERIZATION: A MICHAELIS-MENTEN-LIKE EQUATION 

The average Dwell time can be calculated by substituting (fT"2")) into the definition 

r+oo 

(t}= / tf(t)dt. (17) 
Jo 

Hence, the expression for (t) 

1 1 1 w v ( 1 1 \ Sl„ / 1 1 \ 

* = — + — + — + -7T — + — + -7T + 18 

u)\ u> 2 0J3 oj p + lip \u>4 CJ5 / ujp + S2p \12 4 SZ5/ 
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written in terms of oji,u>2) k^Co^a^, O4 and O5 has a clear physical meaning. In the special case Q p = 0, (i) 



E 



Ul, 



. However, if fl p 7^ 0, then the sum cj 4 1 + oj 5 1 is replaced by the last two terms of ()18[) where w p / (oj p + fl p ) 



and tt p /(ujp + rip) are the weight factors associated with the two paths emanating from the state labelled by 3. 
Finally, in terms of the rate constants for the mechano-chemical transitions in the model depicted in fig|3] 



(t) = i = -(l + ^ 
V uj a \ uj h i 



1 + 



uj r2 



up 
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uip + rip 



1 

UJbf 



1 



Uh2 



1 

UJh2 



1 

OJhl 



UJ r2 



1 

6/ 



1 



Obr 
Oft2 



+ Op 
1 

Oft2 



(19) 



Next we'll show that the equation (|T9|) can be re-expressed in a form that resembles the Michaelis-Menten equation 
for simple enzymatic reactions [28| . Consider an enzymatic reaction of the type 



E 



W +l Win 

■ s£es-¥e- 



p 



(20) 



where -E is the enzyme, 5 is the substrate and P is the product of the reaction catalyzed by E. Given that the total 
initial concentration of the enzyme is [E]o, the rate V of this reaction is given by 



1 

V 

where the maximum possible rate of the reaction is 

'7n aa 

and the Michaelis constant Km is given by 



K 



M 



v m 



v m 



[S] 



W + 2[E\o = 0J2 



K M = 



Ul + l 



(21) 



(22) 



(23) 



In the case of our model, we assume that the "pseudo" first order rate constant uj a can be written as oj a = oj^[tRNA] 
where [tRNA] is the concentration of tRNA molecules in the solution. Treating tRNA molecules as the analogues of 
the substrates in an enzymatic reaction, equation (1191) can be re-expressed as a Michaelis-Menten-like Equation [28j 



<*) = 



1 



K M 



1 



V m 



Vmax [tRNA] 



where 



Vn, 



eff 
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1 



UJ r2 



up -r 

OJbr 
Uh2 



and the Michaelis constant 
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Uh2 
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tt bf 



Qbr 
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eff 
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(24) 



1 



(25) 



(26) 



with 



eff 



UJ rl 



0J r 2 



U> r \UJ r 2 



UJhi UJp + flp 0Jhl{Wp + Op) 



eff 



(27) 



Thus, the mean dwell time for the ribosomes follows a Michaelis-Mentan-like equation. This result is consistent 
with the experimental observations in recent years [67l - [72| that, in spite of the fluctuations of an enzymatic reaction 
catalyzed by a single enzyme molecule, the average rate of the reaction is, most often, given by the Michaelis-Menten 
equation. What makes the Michaelis-Menten-like equation for the ribosome even more interesting is the fact that a 
single ribosome is not a single enzyme molecule, but it provides a platform for the coordinated operation of several 
bio-catalysts. 
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V. EFFECTS OF CROWDING ON THE DWELL TIME DISTRIBUTION 

It is well known that most often a large number of ribosomes simultaneously move on the same mRNA track each 
polymerizing one copy of the same protein. This phenomenon is usually referred to as ribosome traffic because of its 
superficial similarity with vehicular traffic on highways (2(| I29M4H . Suppose, i denotes the number of codons that a 
ribosome can cover simultaneously. Extending the prescription used in our earlier works on ribosome traffic [26j | for 
capturing the steric interactions of the ribosomes, we replace the equations ([5]), (0 and © by 



dPa(i,t) 
dt 

dt 

dPi(i,t) = 
dt 



-- oj bf P 4 {i,t) - u br P 5 (i,t) - u h2 P 5 (i,t)Q(i\i + t) (28) 
Q bf PZ(i,t) - VbrPs (*,*) - n h2 P;(i,t)Q(i\i + £) (29) 



w M fi5(*-M) + n M JE(i-M) Q(yj*-i + (30) 



where Q(i\j) = 1 — P(±\j) is the conditional probability that, given a ribosome in site i, site j is empty. All the 
other equations for P M (z,i) remain unchanged. Note that these equations have been written under the mean-field 
approximation. 

In the limit L — >• 00, all the sites are treated on the same footing so that the site-dependence of P„(i,t) drops out. 
In this limit, Q takes the simple form (2(| 

W + 4 = TT^5- ( 31 ) 

where p is the number density of the ribosomes (i.e., number of ribosomes per unit length of the mRNA track). 

Therefore, in ribosome traffic the distribution of the dwell times of the ribosomes is given by the expression ([12")) 
where uih2 and Qh2 are replaced by L^hiQ and SthiQ, respectively. Using the expression (|3ip. we get the DTD for the 
given number density p of the ribosomes. For the purpose of graphical demostration of the effects of ribosome crowding 
on the DTD, we use the coverage density p cov — pi. In FigJ7]we plot the DTD for a few different values of p cov . The 
higher is the density, the stronger is the hindrance and longer is the dwell time. Moreover, higher coverage density 
introduces stronger correlations which our mean-field equations ignore. Consequently, our analytical predictions, 
based on equations written under the mean-field approximation, deviate more and more from the corresponding 
simulation data as the coverage density increases. 

VI. COMPARISON WITH EXPERIMENTAL DATA 

The distribution of dwell times involves essentially five different rate-determining parameters cj; (i = 1, 2, 3, 4, 5) if 
the translational fidelity is perfect. This is consistent with the earlier observation |24[ that the best fit to the simulation 
data was obtained with five rate-determining parameters. However, the dwell time distribution obtained from the 
single- molecule experiments (more precisely, single-ribosome manipulation) fit reasonably well with a difference of just 
two exponentials [23j . This strongly indicates the possibility that only two of the five rate-determining parameters 
were rate-limiting under the conditions maintained in those experiments. However, for a quantitative testing of the 
predictions of our theoretical model, one should use a mRNA template with a homogeneous codon sequence and only 
two species of aa-tRNA one of which is cognate (corresponding to the codons in the coding sequence of the mRNA) 
and the other is either near-cognate or non-cognate. 

In order to clarify the proposed experimental set up, let us consider the concrete example shown schematically in fig. 
[8j This example is essentially the protocols used by Uemura et al. in their single- molecule studies of translation 
in real time. The actual coding sequence to be translated consists of n c number of identical codons; in figJS] n c = 6 
and each codon is UUU which codes for the amino acid Phenylalanine (abbreviated Phe or F). The coding sequence 
is preceeded by a start codon AUG and is followed by a stop codon UAA. The start codon itself is preceeded by an 
untranslated region (UTR) at the 5'-end of the mRNA; this is required for assembling the ribosome and for stabilizing 
the pre-initiation complex. At the 3'-end, the stop codon is followed by a sequence of n nc non-coding codons UUU 
(n-nc — 4 in fig IB}; this region merely ensures the absence of any "edge effect", i.e., the translation is not affected 
when the ribosome approaches the 3 '-end of the codong sequence. A good choice for the corresponding near-cognate 
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F F F F F F 
■~j AUGUUUUUUUUUUUUT^^ — 3' 

UTR Start Stop 

FIG. 8: A schematic description of a mRNA with homogeneous (poly-U) coding sequence (adapted from refs.[2p|. l2lj|V 

tRNA would be tRNA LeM because it is cognate for the codon CUU which codes for Leucine (abbreviated L) . The two 
distinct species of aa-tRNA should be labelled by fluorescent dye molecules of two different colors which could be, for 
example, Cy5 (red) and Cy3 (green). Equal concentrations of Cy5-labelled ternary complexe aa-tRNA p,le -EF-Tu-GTP 
and Cy3-labelled ternary complex aa-tRNA ieu -EF-Tu-GTP should be made available in the medium surrounding the 
ribosome to avoid any bias arising from difference in their concentrations. The color of the fluorescence pulse identifies 
the amino acid monomer species that elongates the polypeptide by one unit; monitoring the colors of the fluorescence 
pulses one would get an estimate of the translational fidelity <f>. Moreover, the time interval between the arrival of 
the successive aa-tRNA molecules provides an estimate of the dwell times of the ribosome. Thus, using this optical 
technique one would get not only the DTD, but also the fidelity of translation. Usually the coding sequence of such 
poly-U mRNA strands is quite short. Therefore, for collecting enough data to extract the DTD, the experiment has 
to be repeated sufficiently large number of times. 

The effects of kinetic proofreading on the rate of translational error has been investigated experimentally by several 
different methods in the last three decades. However, most of those methods (see, for example, ref. [7^, |73) require 
bulk samples. But, the more recent single-molecule FRET technique used in the study of tRNA selection [ijj seems 
to be more suitable for testing our predictions on the effects of kinetic proofreading on the DTD. 

The cluster of ribosomes translating the same mRNA simultaneously is usually referred to as a ■polysome. For 
studying the effects of ribosome crowding on the DTD, one has to measure simultaneously the DTD and the polysome 
size 17a. 



VII. SUMMARY AND CONCLUSION 



In this paper we have presented a theoretical model of translation that captures all the main steps in the mechano- 
chemical cycle of a ribosome during the elongation stage. This model also accounts for translational fidelity, kinetic 
proofreading and the crowding of the ribosomes on the same mRNA track. In principle, this model can be extended, 
by increasing the number of "chemical" states, to account for some of the sub-steps which have not been treated 
explicitly in the version of the model presented here. 

In spite of the details already incorporated in this model, we have suceeded in carrying out an analytical calculation 
of the distribution of the dwell times of the ribosomes at the successive codons. We have compared this theoretical 
estimate of the DTD with the corresponding numerical data which we have obtained from direct computer simulation 
of the model. If the motion of the ribosome is not hindered by the presence of any other ribosome on the same 
mRNA track, our analytical treatment yields the exact expression for the distribution for the times of its dwell at the 
successive codons. In this case, excellent agreement between theoretical prediction and simulation data is observed 
provided the data are averaged over sufficiently large number of samples. 

However, because of the mean-field approximation made in capturing the effects of crowding of the ribosomes the 
corresponding analytical expression is approximate. Therefore, the higher is the coverage density, the larger is the 
deviation of the analytical estimate from the simulation data which are averaged over many samples. 

We have analyzed the dependence of the DTD on some of the crucially important kinetic parameters of the 
model to elucidate the physical implications of the result. Our results are in good qualitative agreement with the 
experimental data reported in the literature (23l. l24j. However, for the reasons explained in the sections HI and HT1 it 
is not possible to compare these experimental data quantitatively with the analytical expressions of DTD which we 
have reported in this paper. We hope our theoretical predictions will stimulate further experimental studies along the 
lines suggested briefly in section IVII although some technical hurdles may hinder quick progress. A combination of 
the single-ribosome experiments and bulk measurements may be required for comparing the theoretically predicted 
variation of the DTD with the concentrations of tRNA, GTP and other key molecules involved in translation as well 
as with the increase of futile cycles and crowding. 

It would also be desirable to extend our model in future to account for some important features of translation; some 
of these possible extentions are listed below, (i) Capturing the sequence heterogeneity of real mRNA templates will 
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open up larger number of branched pathways in Fig|3j each corresponding to a distinct species (cognate, non-cognate 
or near-cognate) of amino acid. Moreover, at each step of the ribosome, the rate constants for the correct and 
incorrect pathways will also depend on the codon under consideration, (ii) The interaction of the growing polypeptide 
with the exit tunnel in the ribosome and the spontaneous folding of the nascent protein as it comes out of the tunnel 
may affect the rate of translation. It may be possible to capture these effects by making some of the relevant rate 
constants, e.g., ujh2, Qhi, etc. on the instantaneous length of the growing polypeptide which, in turn, is identical to 
the length of the codon sequence already translated by the ribosome. (iii) Maintaining correct reading frame on a 
sequence- homogeneous mRNA (e.g., poly-U) may be more difficult than on an inhomogeneous sequence. It may be 
desirable to extend our model allowing the possibility of frameshift errors (76| although, at present, the prescription 
for this extension is not obvious. The DTD for such a "realistic model" may be obtained numerically because an 
analytical derivation of the general expression may not be possible in the forseeable future. 

Acknowledgements: One of the authors (DC) thanks Joachim Frank and Joseph D. Puglisi for useful correspon- 
dences. 
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Solving the equations ([I])-© under the initial condition (fit))) we get 
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where lji, U2 and W3 are the solutions of the cubic equation (|13p while W4, W5 are the solutions of the quadratic 
equation (fT4| and ^4, ^5 are the solutions of the quadratic equation (fl"5]) . The coefficients (2/ = 1,2,..., 5) and 
D4, D5, which are determined by the initial condition (p~0|) . are as follows: 
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